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ABSTRACT 

We consider the presence and evolution of primordial density perturba- 
tions in a cosmological model based on a simple ansatz which captures - 
by providing a set of effective gravitational field equations - the strength 
of the enhanced quantum loop effects that can arise during inflation. Af- 
ter deriving the general equations that perturbations obey, we concentrate 
on scalar perturbations and show that their evolution is quite different than 
that of conventional inflationary models but still phenomenologically accept- 
able. The main reason for this novel evolution is the presence of an oscillating 
regime after the end of inflation which makes all super-horizon scalar modes 
oscillate. The same reason allows for a natural and very fast reheating mech- 
anism for the universe. 
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1 Introduction 



During the inflationary era infrared gravitons are produced out of tlie vac- 
uum because of the accelerated expansion of spacetime. Such a production 
can only occur for particles that are light compared to the Hubble scale with- 
out classical conformal invariance; gravitons and massless minimally coupled 
scalars are unique in that respect. 

The self-gravitation of the vast ensemble of inflationary gravitons must 
act to slow the expansion rate |1] . This effect is inherently non-local because 
it couples the local graviton energy density with the potential induced by 
the interaction stress of the gravitons throughout the past light-cone. This 
suggests that the relevant effective field equations should be non-local. Non- 
local models of cosmology have been much studied because they can avoid 
the problem that de Sitter must be a solution for any local, stable theory, 
and because non-local couplings between different times can ease fine tuning 
problems [21 [3] . 

Quantum gravitational loop corrections which can be computed during 
inflation grow like the logarithm of the inflationary scale factor [H [5l |6l [7] . 
It should ultimately be possible to derive the most cosmologically significant 
part of the effective field equations by summing the series of leading infrared 
logarithms. Starobinsky has proposed a technique for accomplishing this [S], 
and Starobinsky and Yokoyama have applied it to scalar potential models 
[9]. Starobinsky's method has recently been extended to Yukawa-coupled 
fermions [10] and to scalar quantum electrodynamics [11]. It has not yet 
been extended to quantum gravity but there are reasons for believing that 
some version of it can be [12] . 

Although full control of the effect requires a non-perturbative resumma- 
tion technique, one can attempt to anticipate the results of such a formalism 
in a variety of ways. One approach is to simplify the full quantum gravi- 
tational dynamics by assuming that the exact graviton remains transverse- 
traceless and free, but propagates in the background geometry of an effective 
scale factor which is determined from the expectation value of the Qqq gravi- 
tational constraint equation of motion p^. The simplified theory retains the 
proper perturbative limit for de Sitter spacetime and may provide the basis 
for a tractable non-perturbative formulation. 

Another approach is to use the physical principles responsible for the 
non-trivial quantum gravitational back-reaction on inflation and construct a 
phenomenological model which we can then directly evolve. In a previous pa- 
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per [2] we proposed a phenomenological model which can provide evolution 
beyond perturbation theory. In one sentence, we constructed an effective con- 
served stress-energy tensor Tfj,^[g] which modifies the gravitational equations 
of motion: 

G^!/ = R^iu - ^9fiuR = -^gfiu + SnGT^i^lg] . (1) 

and which, we hope, contains the most cosmologically significant part of the 
full effective quantum gravitational equations. 

Our physical ansatz consists of parametrizing Tf^y[g\ as a "perfect fiuid": 

T^,u[g\ = {p + p)u^Uu + pg^,u , (2) 

with the gravitationally induced pressure given as the following functional of 
the metric tensor: 

p\g]{x) = A^f[-eX]ix) , ^ = 5^ > (3) 
where the function / satisfies: 

f[-eX] = -eX + 0{e') , (4) 
and where the scalar d'Alembertian: 



= d,{g'^''^d,) , (5) 



□ 

' -9 

is defined with retarded boundary conditions. The induced energy density 
p[g] and 4- velocity u^[g] are determined, up to their initial value data, from 
stress-energy conservation: 



which implies: 



D^T^. = , (6) 

u-dp = D^[{p + p)u''] , (7) 

u-dp = -{D-u){p + p) , (8) 

{p + p) u- Du^ = -{dy + UyU - d)p . (9) 



^Hellenic indices take on spacetime values while Latin indices take on space values. 
Our metric tensor 5^1^ has spacelike signature and our curvature tensor equals: R"^^^, = 
r"i./3,;. + r"pp T^^^ -{^^ v). The initial value of the Hubble parameter is iH'^ = A. We 
restrict our analysis to scales M = ( A/87rG')i below the Planck mass Afpi = G^^ so that 
the dimensionless coupling constant e = GA of the theory is small. 
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The 4-velocity is chosen to be timehke and normahzed: 

g^'''u^u, = -l =^ n'^ix^;. = . (10) 

The purpose of this paper is to study the behaviour of the primordial 
scalar perturbations in the phenomenological model summarized above. Due 
to the non-local structure of the model the standard perturbation analysis 
must be extended and this is done in Section 3 after the relevant cosmological 
background is presented in Section 2. The evolution and novel features of 
the scalar perturbations are addressed in Section 4 and their normalization 
in Section 5. A result of the novel behaviour of our scalar mode functions 
after the end of the inflationary era allows a natural and very fast mechanism 
for reheating the universe which is described in Section 6. Our conclusions 
comprise Section 7. 



2 The Cosmological Background 

The large-scale homogeneity and isotropy of the universe selects Friedman- 
Robertson- Walker {FRW) spacetimes as those of primary cosmological in- 
terest; their line element for zero spatial curvature equals in co-moving coor- 
dinates: 

df = g^y{t) (Ix'^dx" = -dt^ + a^{t) rfx ■ rfx . (11) 

Derivatives of the scale factor a{t) give the Hubble parameter H{t) - a mea- 
sure of the cosmic expansion rate - and the deceleration parameter q{t) - a 
measure of the cosmic acceleration: 

(12) 
(13) 

For these spacetimes the stress-energy tensor ([2]) equals: 

Too = uqUq{p + P) - p = p , (14) 
Toi = , (15) 
Tij = UiUj{p + p) + Qijp = pijp . (16) 



m 
lit) 



ciyt) d 
a{t) dt 
a{t) d{t) 



Ina(t) 



HHt) 
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An immediate consequence of isotropy and the normalization condition f llUp 
is: 

The Ricci tensor and Ricci scalar become, respectively: 

3a 



(17) 



and: 



Rij 
R 



, 



d 2a^ 

- + — 

a 



6a 6a^ 

— + — 

a 



3//2 + H\-g^^ 



9ij 



+ 6H 



(18) 
(19) 

(20) 
(21) 



In view of fll4m6[ [T81I2T]) . the gravitational equations of motion ([T]) take 
the form: 



3H^ = A + SttG p 
-2H -3H^ = -A + 8nGp 

while the conservation equation (E]) becomes: 

p = -3H{p + p) . 

The latter implies that: 

1 /■* 



-p{t) = -Pit) + 



dt' a^it') Pit') . 



(22) 
(23) 

(24) 
(25) 



a3(t) Jo 

When acting on functions which only depend on co-moving time, the 
scalar d'Alembertian ([5]) for FRW geometries equals: 



[dl + 3Hdt 



so that its inverse is: 



1 

n 



* 1 
dt' 



0?{t') 



dt" a' it'" 



(26) 
(27) 



Consequently, the source X can be written as follows: 

dt" a%t") \l2H\t") + 6H\t")] . (28) 



X 



1 - /■* 1 

^R = - I dt' 
□ 



/o a3(t') Jo 
Note that we have taken the initial time to be at t = 0. 
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3 Perturbations 



Small deviations from the homogeneity and isotropy of the FRW geometries 
are necessary to address, among other issues, the existence of primordial 
density perturbations. To account for these deviations, we first define small 
but general perturbations of all relevant variables about their FRW values. 
In co-moving coordinates we have: 

(29) 
(30) 
(31) 
(32) 

+ /ioo(t,x) , (33) 
a(t)/ioi(t,x) , (34) 
= a2(t) [5,, + /i,,(t,x)] . (35) 

We then substitute ( I29]|35l) in the gravitational equations of motion ([1]): 





= p{t) + Ap{t,^) , 


p(t,x) 


= p{t) + Ap(t,x) , 


tio(t,x) 


= + Awo(^,x) , 


Ui{t,X.) 


= Ui + a{t) A'Uj(t,x) , 


fl'oo(i,x) 


= 9oo + /ioo(^,x) = -1 




= goi + a{t) /ioi(^,x) = 




= gij{t) + a'^{t) hij{t,x) 



2 ^ 



-A 5^., + SttG 



-A .5^,, + T^,, 



(36) 



and obtain equations to the desired order in the perturbation; for our pur- 
poses, the first order equations will suffice. This is a tedious process which 
to some degree can be simplified by a proper choice of coordinate system and 
field variables. It turns out that conformal coordinates (77, x) and a set of 
gauge invariant variables is the optimal choice [16j. 

In conformal coordinates - which we shall use thereafter in this Section 
- the background invariant element is proportional to that of fiat spacetime: 



dr]^ + dx- dx 



ds^ = g^uiv) dx^dx" = a^{r]) 
The relation between co-moving and conformal times is: 

dt d 
a{t) dr] 



dr] 



d 

a — 
dt 



The corresponding perturbations take the form: 
p(77,x) = p{r]) + Ap(77,x) , 



(37) 



(38) 



(39) 
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p(r/,x) = p{r]) + Ap(r/,x) , (40) 

M^(r7,x) = Uf,{r]) + a{r]) AUf,{r],:)i.) = 0(77) [ -5^° + Au^(77, x) ] , (41) 

9ii,yiV:^) = V(??) + «^('7) = a'^iv)[VtMu + h^^{t,x.)] , (42) 

where we have used the conformal analogue of (1171) : 



u,. 



(43) 



• The Left Hand Side 

As a result of straightforward manipulations, the left hand side of the field 
equations (!36|) can be written as: 



with: E 



a" a' 2 
2 4 — 

a3 



a" a' 2 



-2^ + 



a3 ' 



a" a' 2 
-2- + — 



00 



a" a' 2 
-2- + — 



- ( V"' 5% - 5% v'" ) [ hop,. + ho.,, - h'^„ 

+ \h^^P 4- h ^Pt^ - hP 'P - hP 1 — 
[ " ,pu ^ "'up V p " p V \ 

n 1 

-5^ 



hP'^ _ hP 

,p<T "pa 



1 

2^ ■ 



(44) 



(45) 



(46) 



It is convenient to 3+1 decompose the background value in fH5|) : 

_ _ 

70 



G° 

G^, = , 



a" a' 2 



-2^ + 



a3 



(47) 
(48) 

(49) 



^In conformal coordinates, the prime superscript denotes differentiation with respect 
to conformal time. 
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as well as the first order perturbation /S.G^^ in 

^'2 _ _ ]^ 



-3 —rhQQ + — 



+ 



2a2 



^7 /ij,^; hi-iJ, 



^OO.j + 



2a2 



+ 



a" a' 2 
-2- + — 



^00 + ^ 



/iqq + 2/iofc.fc — /il 



2a2 



00 



a 

+ — 



2^o(ij) + Kj 



(50) 
(51) 



(52) 



:*r Identities 

In addition to the equations presented throughout the main text, various 
of the following expressions have been used to obtain the results of this 
subsection: 



I cO 



a'S 



TP + Ar'' 



',p(7 



n rO rO 



5P AO^ + ^P^^O PO^^^ 



a 
a 



hP + hP - /), 'P 



(53) 
(54) 

(55) 
(56) 



• The Right Hand Side 

We write the gravitationally induced stress-energy tensor T'^ as: 



Its background value T^, is: 

TP = {p + p)u>^u, + pb^^^ , 
and, using fH5]) . its 3+1 decomposition takes the form: 



T° = 



T) = pb,. 



(57) 
(58) 
(59) 



The perturbation AT^^ defined in (1571) equals: 

AT^^ = M^M, A(p + j9) + (p + p) A^M,) + Ap^'t , (60) 
and, in view of tlie 3+1 decomposition is given by: 



AT° 







A^ 



- Ap + (p + - a Am° + a'^ Auq 
(p + p)a"^Awi , 



(61) 
(62) 
(63) 



* The induced pressure deviation Ap 
The starting point for the exphcit computations is the induced pressure 
ansatz ^ which we shall expand to first order about the background ge- 
ometry and determine Ap. Knowledge of Ap and use of the conservation 
equations will allow us to obtain the remaining deviations Ap and Am^, up 
to initial value data, [f] The aforementioned expansion is: 

p = A'f[-eX] = A'{f[-eX] - /'[-eX](eAX) + O(e^) } , (64) 

where: 



Therefore: 



p 
Ap 



A'f[-eX] , X 
-eA^f'[-eX] AX . 



Now the first order perturbation AX equals: 



1 



AX = A{n-^R) = ^ Ai?-(An)X 



(65) 



(66) 
(67) 



(68) 



and we must evaluate: 

(i) the Ricci scalar perturbation AR: 



R = R + AR 



(69) 



■^As we shall see later in this subsection, the perturbation Auo can be computed inde- 
pendent of the conservation equations. 
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for which a straightforward calculation gives: 



R = 6 



AR = 6^ hoo 



a 

1 r 
+ - 



34 



(a) the d'Alembertian perturbation AD: 

n = □ + An 

which directly follows from the definition (jS]): 



-9 



0? ri^"" d. 



dl-2-d. 



V'--doa'do 
a'' 



AC 



A 



+ 



d,. 



+ 



{^9 



-g dv 



d,. 



As a result, the action of AC on the background source X gives: 



(AC) X = ^X'\-h' ~ /i'oo,o + h^,^ 



hoo ^ do a doX 



2 

= ^ ~ 2^' ~ ^00,0 + + -R^OO • 

Subtracting equation (1771) from equation fl7T|) gives AX via fl68l) and 
turn, Aj9 via ([6] 



Ap 



1 



a' r 



eA^ f[-eX]x-^{3^ h'^, - 2ho^, + 



+ 



1 r 



V^hoo - V^hii - 2h'oi i + hf^i + hij^ij 



1 



1 



+ — X' -h' + K 



00,0 "'Oi,i 



ho 



* The deviation Auq 
It is important to note that Auq can be directly determined from the per- 
turbation of the 4- velocity timelike condition ffTOl) : 

g'^'u^u, = -1 =^ {Ag'''')u^u, + 2g^'''u^{Au,) = (79) 

^ -/loo + - Alio = . (80) 
a 

We trivially conclude that: 

Auo = ^ahoo , Am° = ^ ^oo • (81) 
As a result, AT°q - given by equation fl6T|) - is simplified: 

Ar°o = -A/9 . (82) 



:*r Identities 

In addition to the equations presented throughout the main text, various 
of the following expressions have been used to obtain the results of this 
subsection: 



h 



-h 



00 



hi. 



Av^ 



g = a 

h 







1 


1 


V-9 





- h , 
2 

5\X' . 



Ag 
A 



2a^ 



h 



When acting on functions that only depend on conformal time: 

1 . 1 111. 



do a? do 



(83) 
(84) 
(85) 
(86) 
(87) 



(89) 



• The Conservation Equations 

Of the conservation equations ([3inD, we shall use ([7]) and to express Ap 
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and Am„ in terms of Ap. As a first step, we write them more explicitly: 



-g {p + p) 
(p + p) u^.^ u 

and study them up to first order. 



(90) 
(91) 



(i) Zeroth order 

Because the background 4- velocity is given by fH3l) . the background value of 
(I90|) equals: 

(92) 



a^p' 



and integrates to: 



1 n 



p + p = I dr] a? p 



(93) 



which is the conformal analogue of equation (!25l) . The background value of 
(I9T1) leads to a tautology of the form = 0. 



(ii) First order 

When the perturbations defined in fl39H42p are applied to the conservation 
equations (190119 ip they respectively lead to: 



(Av^^) {p + p)u^ + [a{p + p)) + (p + p) {A 
= {A^)u^d^p + {Au^')d^p + v^M^9^(Ap) , 



(94) 



(A(p + p)) u^.^^ + {-p + p) {Au^.^J) + {-p + p) u^.^, {Au^) (95) 
{Ag^''') + {Au") + {Au^)v!'] d,p - \g^''' + u^v!'] d,{Ap) . 



Use of various background values and identities reduces equations fl94tl95p as 
follows: 



^0 a^A{p + p) 



a\p + p)[di Au' + ^KA + Ap' , (96) 



a\-p + p)UAu^y + 2-Au^^--5%Au^ + - ^ VJ } (97) 

= a' p' - [ Au^ + Wp - a' (?7^" - 5\) d,{Ap) . 
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The last equation can be 3+1 decomposed and further reduced: 



a%p + p) (a^ Au') 



(98) 



a^{p + p)hoi\ + - a^{p + p) dihoo 



-a^d,{Ap) , (99) 

where we have extensively used (!92|) . Equation (!98|) is trivially satisfied by 
the solution ( jSTl) for Au^. Given the perturbation Ap, equation (!99l) deter- 
mines Au* up to its arbitrary initial value. Then, equation ( j96ll has enough 
information at its disposal to determine Ap, again up to its arbitrary initial 
value. Therefore, AT'^ has been completely specified to first order. 



Identities 

In addition to the equations presented throughout the main text, various 
of the following expressions have been used to obtain the results of this 
subsection: 

(100) 



2 " f 1/ ) 



Auf, 
Aii^ 

d,^P 



up 



(Agf^nuu + r" (Am,) 



Au^^^^ + T% (AtiO + (Ar 



up. 



a{Au^ 



5\P' 



6" — 



«^(P + P) 



, (101) 
(102) 
(103) 
(104) 

(105) 
(106) 



• The Full Equations 

At the cost of being redundant, we list the results derived in this section for 
the gravitational equations of motion (l36l) up to first order in the perturba- 
tions (12911421). 

(i) Background equations: We are given a background FRW spacetime char- 
acterized by the scale factor a, and an induced background pressure p ; then: 

.'2 



(00) 



-3^ 



-ml 



SttG 



1 



p -I drj p' 



(107) 



12 



= 



a" a''^ 
-2 — + — 



a** 



S.. 



(108) 

'iHldij + 87rGp% . (109) 



(a) First order equations: We are given the spacetime perturbation h^y and 
the induced pressure perturbation Ap ; then: □ 



(00) 
(Oz) 



AG°o = 87rG'(-Ap) 



SttG (p + p)- Aui 
a 



(110) 
(111) 
(112) 



• Scalar Perturbations 

The very nature of the gravitational equations of motion allows for scalar, 
vector and tensor perturbations. Of these, it is the scalar perturbations that 
have, up to now, the biggest phenomenological interest. A general scalar 
perturbation is characterized by four scalar functions (p, B, i/j, E which are 
defined as follows [T6l: 



/loo = -20 , 
h^i = —Bi , 
hij = -2ifj6ij 



2E 



(113) 
(114) 
(115) 



and which can be combined into two gauge invariant scalar field variables 



a{B- E') 



-{B- E') 
a 



'B' - E") 



^ ^ ^ + ^m-E') 
a 



(116) 
(117) 



To deduce the equations that are obeyed by the scalar perturbations we 
must isolate their contribution to the full gravitational equations using (11131 
11151) and then express it in terms of the invariant variables (11 16111 17|1 : 



^To avoid prohibitively long expressions we have not expressed all variables in terms 
of hfj^i, and Ap. 
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(i) The left hand side 
The 3+1 decomposition of {AG^j^)s equals: 

2 _o r , a' 



AG* 



Us 



— V 



a' 
6^ 



a 



i^di 



^l) + -iB-E' 
a 



d,. 



+ - 

a 
a 



+ (G\-\G\\d,{B-E'^ 



a" a' 2 



+ 



2^' + 0'l + ^" 



(118) 
(119) 
(120) 
(121) 



+ _V2[0 - V; - (fi' -£;")- 2 -^(S - ^') 



+ - ^ - + 5' - E" + 2-{B- E") 



(122) 



J v'($ - ^) + ^" + - (2^' + $') -\ 
[ 2 a 

X 5,,- + 1 (^,,, - $,,,) + &'AB - F) . 



a" a' 2 
2 ^ 



(123) 



(^wj The right hand side 
A similar analysis must be made to derive (AT'^)^. We need only compute 
the scalar perturbations contribution to Ap since the conservation equations 
will give the similar contribution to Ap and Atij. [f] Instead of calculating 
directly from (ITHj) . we shall perform the computation sequentially by first 
considering the scalar perturbations contribution to (!7T|) : 



a" 

-12-, 



b' + ?,ip' -V'^{B - E') 



V2(2V; - 0) - 3V^" + V2(5' - E") 



(124) 



V2(2^-<l>) 



6 r 



2R4> + — {a^ R)' {B - E') + 2R{B' - E" 



(125) 



^The perturbation Amq does not enter {ATIf)s and, at any rate, has been calculated: 
Auq — —a(j). 
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and the similar part of ( 17^ : 



(An)x 



-2R(i) + ^X' \?,iP' + -V^{B-E') 



(126) 



-2R(t) + —X' + - 



-J2 



(B - E') 



~{B' - E") + {B" - E'") -V\B- E') I . 



(127) 



By exploiting ( !73l) . the last term in (I127P can be rewritten as follows: 

1 _ 



[B - E') 



□ 
□ 



X' (B - E') 



X'{B-E') 



+ \(do + 2-)do 



X' (B - E'] 



;i28) 



X' {B -E') - 4 («' R)' {B - E') - 2R {B' - E") 



+ X'\ -iB" 



E" 



2^(B'- E") 



a" a' 2 
2 2 — 

,3 ^ ,A 



[B - E') 



where in the second step we used another direct consequence of fl75|) : 



(129) 



1 



^,do + 2-)doX 



R 



X'' 



a 



2-X' - a^R , (130) 



a 

^ a ^ a 

Appropriately combining equations (11251 11271 11291) . we get the scalar pertur- 
bations part of the source deviation: 
1 

/\R- (AD) A 

s 



(ax 



s n 

= X'{B- E') + 



AR-{An)X 

1 f 2 
□ I a 



2 V2(2^ - $) 2 



6 r 
a 



a 



^" + - (3^' + $') 



-^X' (3^' + $') 



(131) 



The induced pressure first order scalar deviation follows from ([6 

1 f 2 



Aj9 = p'{B-E') - eA^ f'[-eX] x 



, V\2^-<I>) 



6 



a 



^X'(3^' + $') I . (132) 
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(in) The conservation equations 
The scalar perturbations part of the conservation equations 
comes respectively: 

a^A{p + p)]' = i^Ap' + {p + p)[3'ilj' -V^[^Au + B - E')]^ 



be- 



a%p + p) 



Au^ 



a[p + p) -<p 



O' (P + P) 



a\p + p)B 



Ap+ {p + p) 



(133) 

(134) 
(135) 



As before, equations 0133^ I135p determine Au and Ap while equation fll34p 
is automaticaly satisfied by Auq = —acf). 

(iv) An observation 
Consider the scalar perturbations part of the (ij) first order equation of 
motion: 

(AG*^.)^ = 8nG{AT))^ = SttG (Ap) 5V . (136) 

When i ^ j, equation fll36p simplifies considerably, since all terms propor- 
tional to vanish, and allows us to conclude that $ 



^ j] 







^: 



(137) 



where we have used equations fll23[ I132p and the conformal time derivative 
of the [ij) background equation of motion: 



SttG T 



StiG p' 5' 



(138) 



Identities 

In addition to the equations presented throughout the main text, various 

of the following expressions have been used to obtain the results of this 
subsection: 

h = -hoo + ha = 2 4> - Qi) - 2V^E , (139) 

ha = -6^ - 2V^E , (140) 

ho,^, = -W^B , (141) 

/i,,,,- = -2^,i - 2V^E^i , (142) 



h 



-2VV - 2 V^S , 



(143) 



^We have decomposed the perturbation Aui into its transverse and longitudinal parts: 
Ui = Auf + di Au. The transverse perturbation is of no concern to us. 



16 









" a" a'2 " 
a 




-E') - 


a' 
a 


(B'- 


E") - {B" 






" a" a'2 " 

a 


{B- 


o! 

-E') + - 
a 


{B'- 


E") , 
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a? 
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E') , 



(145) 



= E") -2-(B-E') , 

CL 

r n" 

3^' + $' = 3ip' + (f>' + 

-{B"-E"') , 



a" a' 2 
2--2 — 
a 



2(?' 

{B -E')+ — {B' - E") 
a 



a a 



6 



(o^ (5 - E') 



-J2 



G\ - - = — 



a 



-,/// 



a' a" 



,/3 



+ 2- 



a 



3 ' 



^1 „ 0,"' ^ „ a' a" 
i?' = 6 — - 18 —— 



a" 



a'R)' = 6 — 



a' a" 



P 



-eK'X'f[-eX] , 



a 



(146) 
(147) 
(148) 

(149) 

(150) 

(151) 

(152) 

(153) 

(154) 
(155) 

(156) 



• The Invariant Completions 

There is an elementary way to redefine the Einstein tensor perturbation 
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AC^^ and the stress-energy tensor perturbation AT'^ so that both sides of 
the equations of motion are gauge invariant 



AG^^ = STrGATt A^^^ = SttGAT^ . (157) 

The idea is for the differences between (AG^^ , AT^^) and (A^^^ , AT^) to 
simultaneously make the latter gauge invariant and obey the gravitational 
equations of motion. This can be accomplished in general with the redefini- 
tions: 

Ag% = AG\ - GX{B-E') , (158) 

1 



A^° = AG', - (G% - - G^^) {B - E'), , (159) 

AG'. = AG' J - G'.{B-E') . (160) 

For instance, in the case of scalar perturbations these difference terms are 
precisely the only terms of (AG^)5 that are not gauge invariant as can be 
seen from flll9[ I12H I123p . The corresponding redefinitions of the stress- 
energy perturbation are obvious: 

AT^o = Ar°o - f'^\B-E') , (161) 

AT? = Ar° - (n - \ {B - E'), , (162) 

AT) = AT) - T-{B-E') . (163) 

(i) The scalar perturbations equations of motion 
Returning to the scalar perturbations that concern us, we have: 

(00)5 =^ Ag\ = ^'KG{AT\-f\'{B-E')] (164) 

=^ -^V'$ + 64'^'' + = 87rG(-Af) , (165) 

(0^5 =^ AG\ = 87rG{AT° - (n " ^ n) " ^0.} W 

=^ - $'-2^$ = 8^Ga-^(AW) , (167) 

{tj)s =^ Ag'^ = 87^G{AT) -f)'{B-E')} (168) 

^ ^<l>" + 64'^''+ [4^-2^1 $ = 87rG(AP) . (169) 



^Since we are only concerned with scalar perturbations, from now on we shall not 
indicate this explicitly. 
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In arriving at equations ( 11651 11671 1169p - besides the equality $ = \1/ - we 
used expressions (II 19^ I12H I123p respectively, as well as the definitions: 

/\£ = /\p - p\B - E') , (170) 
AW = An + a{B-E') , (171) 
AV = Ap - p' {B - E') . (172) 

Our induced pressure ansatz provides AV: 

-1(9 6 r An' id- 1 

AV = -eAVheXlxi 4V'$-^ -^X'$' (173) 



a 



The remaining two quantities A£, AlA are obtained from the conservation 
equations ( 11751 11761) . 

(a) The initial value problem 
The appropriate set of initial value data for the system of equations of motion 
( 11651 1167[ I169P consists of the following: 

{A8\=^^ & {AU\=^^ : unrestricted , {AV\=^^ = . (174) 

The requirement {AV)no = comes from the perfect fluid form ([2]) of the 
stress-energy tensor and the induced pressure ansatz (jS]) ; is the retarded 
Green's function and vanishes on the initial value surface. The initial value 
data {A£)ng, (AW)^^ are free and via equations (I165[ I167P determine ^r/o^ 
Then, equation (I169P determines 

(Hi) The scalar perturbations conservation equations 



The equations of motion are augmented by the conservation equations (I133[ 
I135p whose invariant completions are: 



AS 



a^p + p)AU = -a"" AV - a^ip + p)^ . (176) 



[a 



AU. 



AV + 3a3 (p + p) - (p + p) ) , (175) 



(iv) The dynamical content 
The fact that (Af)^^ & {AU)^g are unrestricted and, therefore, "^"^^o 

^We have again decomposed the perturbation AUi into its transverse and longitudinal 
parts: AU^ = AUf + 9, AU. 
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also unrestricted, implies that there is a scalar degree of freedom which be- 
comes dynamical due to the presence of our gravitationally induced stress- 
energy tensor T^j^[g\. 

This is to be contrasted with the situation where the origin of the stress- 
energy tensor is the matter sector of the theory; if we denote the relevant 
deviation of the latter tensor by AG^^,, we have: 



(00)5 



2 „9 ^ o' r I , a' ^ 



87rGAe°o , (177) 



As. 



a 



SnGAQ^^ = SivGd.iAQ) , (178) 



where the last equality in fll78p is true because we consider scalar perturba- 
tions. From fll78p we immediately deduce: 



$' H $ 

a 



— StvGAQ . 
2 



Substituting ffTT^ in the (00) ^ equation f lTTTD : 
2 



S-SirGAQ 
a 



SnGAQ^ 



) 



(179) 



(180) 



allows us to conclude that $ has no dynamics since we can solve for it from 

(HHoD. 

Identities 

In addition to the equations presented throughout the main text, various 
of the following expressions have been used to obtain the results of this 
subsection: 



-L < 



rpl 



p'5' 



-S-iP + p) , 
a 

° 3 ^ 



(P + P) 



(181) 
(182) 



4 Scalar Perturbations Equation Solutions 

In this section we will investigate solutions of the non-local evolution equa- 
tion fll69p . It is most convenient, for this purpose, to return to co-moving 
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coordinates and we shall do so: 



2 f'[-^X] 

■ ur X — - X 



f 

J cr 



1 \ f 2 ^ 

W $ + 5if$ + -X$ + — $ 
n/ 1 3 3a2 



where the critical point X^r and oscillation frequency u) are 

1 - 87re/[-eX] = , a;^ = 24^e2A/;^ 
The scalar d'Alembertian acting on a general function equals: 

= + 3i/9t + ^ . 



;i83) 



(184) 



(185) 



It should be clear that equation f ll83p cannot be solved exactly and we need 
to develop a methodology that will allow us to extract the time evolution of 
the scalar perturbations. 

• Strategy 

We can realize the inverse d'Alembertian using the mode functions u{t, k) and 
u*(t, k) of the massless, minimally coupled scalar which obey the equations: 



u{t,k) + 3H{t)u{t,k) + 



k' 



u(t, k) 







u{t,k) u*{t,k) — u{t,k) u*{t,k) = ia ^(t) 



(186) 
(187) 



The solutions for general a{t) are quite complicated [2U] but we shall only 
require the asymptotic forms long before and long after first horizon crossing: 



k > H{t)a{t) 
k < H{t)a{t) 



u{t, k) 
u{t, k) 



1 



'2k a{t) 
H{h) \ - 
V2F 1 - 

i k^ 

where tk is the time of first horizon crossing: 

k = H{tk) a{tk) 



exp 



+ik 



+~ dt' 



dt' 



a{t') 



Oik^ 



188) 



;i89) 



(190) 
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THE FIVE REGIMES 



before 1st after 1st oscillations matter radiation 



Figure 1; The regimes during which the evolution of scalar perturbations was studied. 



We can construct the retarded Green's function of □ using the mode 
functions u{t,k) and u*{t,k). Hence the action of — D"^ on some function 
fit, k) gives: 



i dt' a%t') u{t,k) u*{t',k) -u*{t,k) u{t',k) f{t',k) ,(191) 

Jo L J 

and further action on — □/(t, k) only gives back the function /(t, k) up to 
homogeneous solutions: 



-g){-n/}(t,A;) = f{t,k) + au{t,k) + Pu*{t,k) . 



(192) 



The constants a and /3 can be expressed in terms of the initial values of / 
and its first derivatives: 



a 



f3 



u*{0,k) f{0,k) -u*{0,k) f{0,k) 
u{0,k) f{0,k) - u{0,k) f{0,k) 



(193) 
(194) 



Our approach will be to divide time evolution into successive regimes and 
obtain reliable approximate solutions to 01831) within each regime. The choice 



of these regimes is dictated by the actual physical evolution of the system and 
by our desire to reliably approximate (11831) . For instance, the last term can 
be irrelevant or important depending on whether the particular mode with 
wave number k has or has not experienced first horizon crossing respectively. 
We shall assume that the wave number k lies in the range for which the mode 
experiences first horizon crossing during inflation, but close enough to the 
end of inflation that the mode is at a cosmologically observable scale today. 
There are five epochs during which we seek approximate solutions to (11831) . 

The strategy is always the same. We first determine the "/' ratio" in the 
left hand side of f ll83p . then we approximate X{t) in terms of H{t). We next 
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identify the "large" part of the curly bracketed term in the right hand side of 
(11831) and extract it from the inverse d'Alembertian using (11921) . along with 
the appropriate homogeneous solution. 



• The Inflationary Regime Before First Horizon Crossing 

The simplest epoch is the first one during which: 

and where the "/' ratio" is very much less than one for most choices of the 
function f{x). The "large" part of the curly bracketed term is so we 

extract it using f lT92l [T85|) : 

-^){4 + 5ff* + ?A-* + |l*)(,*)« 

A;) + {homogeneous) + (^-g) | + ^'^'^|(^'^) (196) 

where (homogeneous) stands for a linear combination of the homogeneous 
solutions u{t, k) and u*{t, k) given by f llSSp . Because of conditions (11951) and 
the form of (11881) . the entire right-hand side of (11831) is negligible and the 
equation reduces to: 

^ + AH^ + (3i/2 + 2H)^ ^ Q . (197) 
Approximate solutions of (I197P are: 

$i(t,/c) ^ ^ , ^v{t,k) ^ , (198) 

a[t) a'^[t) 

and the recognition of this fact allows us also to estimate the first corrections 
from the non-local term. Suppose, for instance, we are correcting the $i 
solution, in which case the non-local term is approximately given by the 
retarded Green's function (11911) - constructed out of (11881) - acting on the 
dominant term ^'^'i: 

(_^)|i + 5//i + |i* + |l*|(a) 
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1 /■* 

—TT / dt' a^it') sin 
ait) Jo ^ ' 



k a{t) 
1 



dt' 



k 



3a{t) Jo a{t') 
1 



3a{t) 



cos 



k 




3a2(f) a{t') 



sm 



(199) 
(200) 
(201) 



The leading correction to $i derives from the non-oscillating term ^ in 
(12011) . However, this correction - which grows like u}'^H~^ In a x $1 - is of 
the same kind as the correction to $1 we would get if we take into account 
the 2if$ term in equation (11971) . |^ Moreover, the dominant correction from 
the oscillating term in (120 ip is, again, computed by inserting it as the source 
in the evolution equation (I183p and noting that the leading time dependence 
comes from the oscillating source term itself: 



A^tk) = Ut,k) X LL-^ x(--)^ COS 



* dt' 

a(f) 



(202) 

The correction is very small long before first horizon crossing, owing to the 
factor of k~'^uj'^a?{t). However, we see that the oscillatory term does grow. 



• The Inflationary Regime After First Horizon Crossing 

During this regime we are close enough to the end of inflation that, for 
instance, the "/' ratio" is nearly unity: 

^ « H\t) , uj^ « H\t) , X{t) ^ -mt) , 4fr^ " 1 ^203) 
a\t) j'\-eX^r\ 

This is a difficult epoch to understand from first principles so we had re- 
course to explicit numerical studies. These revealed an end to the fall off in 
$(t, A;) which characterizes the previous epoch. In fact, the solution changes 
sign and its magnitude seems then to grow slowly. One can understand why 
this happens from the constant homogeneous solution which is built up by 
integrating the retarded Green's function over times t' long before first hori- 
zon crossing. For this case we approximate u{t^ k) by expression (I189p and 
u{t',k) by expression (I188p . We further assume ^{t',k) ^ a^^{t'), and that 

^During inflation H « — 
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the ^^{t') term dominates: 



1 



dt' a?{t') sin k 
(if 



3a? ) 
~ df ' 



cos 



k 



a{t') 



— cos 



a(t") 
k 



3a2(f) a(t') 
^ df " 



a{t') 



C2{k) 



(204) 
(205) 



and we see that the constant is positive: C2{k) > 0. Thus the approximate 
evolution equation during this epoch is: 



-uj^C^ik) . 



(206) 



Even though u"^ ^ the term on the right is not zero, so the a~^{t) fall 

off of ^{t,k) cannot persist indefinitely. When this rapid time evolution of 
$(t, A;) comes to an end the time derivative terms become insignificant and, 
because —H{t) <^ H'^{t), we have: 



3i/2(t)$(t,A;) ^ -u^C2{k) 



$(t,A;) 



uj^C2{k) 



(207) 



Hence the solution changes sign and, because H{t) decreases slowly during 
inflation, the magnitude of the solution increases slowly. That is exactly 
what the numerical simulations show. 



• The Oscillatory Regime 

During the epoch of oscillations the "/' ratio" is still unity, and we also have: 

^ « H'{t) « \H\ « , X{t) ^ -6H{t) , 4!r^ ^ 1 (208) 
a^{t) f'[-eXcr\ 

Expression fl207p implies that $(t, A;) must begin evolving again at the end 
of inflation, so that its time derivatives are no long negligible. This has 
two consequences: first, the non-local term receives substantial contributions 

^'^The contribution to the integration in (|204p from to t is subleading relative to that 
from to tk- 
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from times t' after criticality; and second, the "large" term is $. Therefore 
we can write: 

+ ^^]it,k) ^ ^t,k) + Csik) , (209) 




3" * ' 3a2 J 

where Csi^k) consists of C2{k) plus a new, time independent contribution. 
Employing fl208p . the mode equation becomes effectively: 

^{t, k) + m{t) ^{t, k) + u"^ ^{t, k) ^ -u^ C^ik) (210) 

^ Htk) « -cm + ^'^w^'-'i^^^l+fewi , (211) 

where At = t ~ tcr and $3 (A;), 03 (fc) are constants. Had we included the 
first correction to the right hand side of fl209p the result would have been 
to slightly modify the rate of fall off in the oscillatory term, but it would 
not change the frequency. We emphasize that all the super-horizon modes 
oscillate at this same frequency. This is a profound distinction between the 
scalar perturbations of our model and those of scalar-driven inflation, and it 
has important consequences for reheating. 

• The Subsequent Matter and Radiation Domination Regimes 

Since the field $ couples to matter universally with gravitational strength, 
its oscillations with frequency u will most likely excite the particles with 
masses m u. These will be very heavy particles which will behave - even 
after their excitation - like non-relativistic matter. The very heavy particles 
that were the primary receptors of the energy from the oscillating field $ will 
quickly decay into a hot radiation dominated universe. 

The analysis for these two epochs of matter is much the same. The only 
difference concerns the approximation we make for X{t): 

matter =^ X{t) ^ -2H{t) , (212) 

radiation =^ ^{t) ~ . (213) 

These differences shows up only in the correction to (12091) . which affects the 
rate of fall off but not the oscillatory frequency. Furthermore, more evolution 
also affects the constant parts accumulated from the homogeneous solution: 

^/ ,N ^/,N ^Jk) sinlu At + (bJk)] 

matter =^ $ t, A; ^ -CJk) + ^— Ti^^ ^ 214 

a^(t) 

, ,N ^rJ k) sinl CO At + (br,(k)] 

radiation =^ $ t, A; ^ -C^ik) + ^ , . 215 

a^[t) 
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We have not included the very significant dechne in amplitude which must 
occur due to the flow of energy from the scalar modes into normal matter. 
It seems clear that this will continue until the amplitude of oscillation is 
driven to nearly zero. The final signal for the power spectrum resides in the 
constant C^lk) whose normalization we cannot fix in the absence of canonical 
quantization. 

5 The Normalization of Perturbations 

Our effective field equations govern the time dependence of perturbations. 
We cannot actually solve these equations exactly but let us suppose, for the 
purposes of this discussion, that we could. For each wave vector k that would 
determine two linearly independent solutions, k) and $2(^, k). The full 
content of the effective field equations has been exhausted by expressing the 
perturbation operator k) as a linear combination of these two solutions: 

$(t,k) = ai(k) X $i(t,k) + a2(k) X $2(t,k) . (216) 

We can say what the operator coefficients, Q!i(k) and Q;2(k), are in terms of 
the initial values of ^(t, k) and its first derivative, but the field equations 
alone do not define how these operators commute. That information would 
ordinarily derive from applying canonical quantization to a Lagrangian, but 
in our non-local cosmological model we have no Lagrangian. We must in- 
stead regard the missing information as a separate assumption, which can be 
specified however we wish. It is a fundamental part of the definition of the 
model, every bit as much as the effective field equations were. 

Before stating this assumption, let us clarify the issues in the very simple 
context of a 1-dimensional point particle whose position q{t) obeys the simple 
harmonic oscillator equation: 

q{t) + oj^q{t) = . (217) 

This is a trivial equation to solve, and we can use it to express q{t) in terms 
of its initial values go ^-nd %: 

Q{t) = Qo cos{ojt) + — sin(co't) . (218) 

OJ 

By decomposing the oscillatory functions into positive and negative frequen- 
cies we can identify linear combinations of the initial value operators which 
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must lower and raise the energy: 

9W = o(?o + — je + -(go 



"i/,go±^l = T^cUi^). (219) 

Relation (121 9p is as far as one can go using only the equation of motion fl217p . 
We do not know how go and go commute, nor do we know how the Hamilto- 
nian depends on them. Indeed, these two issues are intimately related. If we 
ignore possible operator ordering ambiguities, the two Hamiltonian evolution 
equations: 

% % dH 

go = - [if, go] = -^^bo,go] , (220) 

are consistent with any solution of the form: 

H = F{8) and [go, go] = , where 8 = ^g^ + ^uj^^^ ■ (222) 

The equation of motion (\217\\ cannot tell us what the function F{S) is. 

Note that there is still an ambiguity even if we assume F{S) is linear - 
which assumption might seem reasonable (athough not necessary) in view 
of the fact that the equation of motion is linear. The ambiguity rests with 
the proportionality constant: any function of the form F{S) = KS would 
reproduce the canonical operator equations fl220ti22ip . Therefore, if we 
write g(t) as a linear combination of canonically normalized creation and 
annihilation operators, the amplitude with which they appear contains a 
factor of the arbitrary constant K: 



q{t) = ^^^{ae— * + aU^^'} , [a, a^] = 1 . (223) 

It is canonical quantization of the simple harmonic oscillator Lagrangian 
which would ordinarily fix this constant: 

L = -mq'^ — -mu'^q'^ =^ K = m . (224) 



There is no operator ordering ambiguity for any linear ansatz. 



28 



In our case, the effective field equations are not associated with any La- 
grangian. Indeed, the simultaneously causal and non-local nature of our 
equations precludes their derivation from any single-field Lagrangian. El A 
correct derivation from fundamental theory would involve an effetive action 
of the Schwinger-Keldysh type in which more than one quantity stands for 
what will eventually be the same dynamical variable, and one obtains the 
equation of motion by varying first and only then setting the different quan- 
tities equal [21]. 

In the absence of an action principle we are forced to regard specification 
of the constant K as an independent assumption, with the same status as the 
effective field equations themselves. This assumption can be made however 
we wish as part of how we define the model. In this context it should be 
noted that no principle seems to preclude the constant depending upon the 
magnitude of the wave vector, k = ||k||. Efl It is immediately obvious that 
we can enforce both the observed magnitude of scalar perturbations, and 
their approximate scale invariance, if only the super-horizon mode functions 
freeze in to constants in the late time regime. We saw in the previous section 
that the super-horizon mode functions do become constant during radiation 
domination, so that is the non-trivial check of our class of models, not either 
the magnitude nor the approximate scale invariance of scalar perturbations. 

These comments need not apply to tensor perturbations. The equations 
which describe the latter are the same as those of general relativity, pro- 
vided the expansion history a{t) is fixed. We can therefore invoke the usual 
canonical normalization of the graviton creation and annihilation operators 
even though that is not, strictly speaking, required. Doing so makes the full 
tensor power spectrum a prediction of our model [17]. 

6 Reheating 

We have seen from Section 4 that the evolution of scalar perturbations reveals 
a profound difference between our model and any scalar-driven model of 
inflation: 

^^Although a partial integration "trick" can give causal non-local field equations, the 
resulting models inevitably suffer an undesirable renormalization of the effective Newton's 
constant which the present class of models avoids by construction [2] . 

^■^Invariance under spatial rotations rules out more general dependence. 
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• In scalar-driven inflation the scalar mode functions oscillate and decay 
until horizon crossing, after which they approach constants; whereas 

• The scalar mode functions of our model decay until horizon crossing, 
then they become approximately constant until the end of inflation, 
after which they oscillate - all with about the same frequency u - until 
enough energy has been dumped into the matter sector to support a 
radiation dominated cosmology. 

Another key distinction is that, whereas no one knows or cares how long 
scalar-driven inflation persists beyond the 50-60 e-foldings needed to solve 
the horizon and flatness problems, our model requires the vast number of 
(GA)^^ > 10^ e-foldings of inflation. This implies that the number density 
n of super-horizon modes at the end of inflation is staggering: 

^ ~ S '''P(^) - • ^^^^^ 

Now consider what must happen when all of these modes start oscillating. 
There is not much energy in any one mode, and they couple only gravitation- 
ally to ordinary matter, so energy flow from them is very weak. However, 
the number density (12251) of modes which participate is so enormous that 
reheating must be practically instantaneous. 

It is interesting to contrast reheating in our model with the way it works 
in traditional single-scalar inflation. In those models the single inflaton zero 
mode undergoes oscillations at the end of inflation. At this stage it becomes 
necessary to assume that the inflaton couples to ordinary matter because 
having the its kinetic energy transferred to gravity and, then, from gravity 
to ordinary matter, would be far too slow to reheat the universe. What would 
happen instead is that the energy would be reshifted away before it had a 
chance to accumulate and thermalize. 

Direct couplings between the inflaton and ordinary matter can result in 
significant reheating. However, it comes at the price of more fine tuning 
because such couphngs cause loops of matter quanta to induce Coleman- 
Weinberg terms in the inflaton effective potential Ks- These contributions 
would not be in terms of naturally small parameters, they would necessarily 
involve the relatively large coupling constants of ordinary matter. For exam- 
ple, a cubic coupling between the inflaton ip and either a spectator scalar x 
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or a fermion field -0 would induce |18j : 



- 9x "PX^ 




(226) 



(227) 



for some renormalization scale jjL. Either of these contributions would render 
the inflaton potential far too steep for successful inflation; the fermionic 
contribution f l227p would actually make the inflaton unstable. Hence these 
effects must be canceled by adding carefully chosen terms to the classical 
action. As long as one is not restricted to renormalizable inflaton potentials 
this can be done, but it represents a new level of fine tuning. 

This additional challenge for scalar-driven inflation derives from the fact 
that reheating draws its energy from the oscillations of just a single scalar 
zero mode. Thus, the inflaton must be directly coupled to ordinary matter 
to give efficient reheating. By contrast, the reheating in our model draws its 
energy from the vast reservoir of super-horizon modes which are naturally 
accumulated during the long epoch of inflation. Because so many modes 
participate, it is not necessary (or even possible) to introduce a new, direct 
coupling to ordinary matter; gravitational couplings will suffice. 



The phenomenological model considered in this study is solely based on the 
graviton and therefore, one would argue, it should have major problems in 
reproducing any realistic scalar density perturbations; after all, the graviton 
is a tensor field possessing four unconstrained initial value data which result 
in its two physical polarizations. Unless we are willing to invoke graviton 
bound states, there is simply no physical scalar degree of freedom present. 
However, this argument ignores the presence of the gravitationally induced 
non-local source term in the field equations. As we showed in Section 3, its 
presence changes the dynamical content of the theory and, besides the two 
graviton polarizations, the scalar $ emerges as a physical degree of freedom 
possessing two unconstrained initial value data: A£^o and AWq . 

Deriving the evolution equation for the scalar $ is a non-trivial exercise 
described in detail in Section 3. The (approximate) solutions to the equation 
in five successive regimes of evolution from the onset of inflaton until late 



7 Epilogue 
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SCALAR MODE FUNCTIONS BEHAVIOUR 


EPOCH 


STANDARD 


MODEL 


from: onset of inflation 
until: 1st crossing 
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falling 


from: 1st crossing 
until: end of inflation 


constant 


constant 


oscillations 




oscillating 


matter 


constant 


oscillating 


radiation 


constant 


constant 



Figure 2; The mode functions evolution in scalar-driven vs gravity-driven cosmology. 

times, sliowed distinctive differences with the standard inffationary picture; 
they are qualitatively recorded in Figure 2. What counts is agreement with 
measurements and what is measured is correlations between different portions 
of the sky. These correlations are non-zero for the scalar part of the generic 
perturbations in the gravitational theory. Scalar-driven inflation makes a 
non-zero contribution dictated by the form of the inflaton Lagrangian. Our 
gravity-driven model also makes a non-zero contribution which - with the 
proper normalization choice - is consistent with the observed magnitude and 
approximate scale invariance of the scalar spectrum. 

The novel feature of the gravity-driven model is the presence of an era 
subsequent to inflation during which a// modes of $ oscillate. By all modes we 
mean all infrared modes since $ is a dynamical degree of freedom emerging 
in the infrared sector of the theory. Now the "receiver" of the energy gener- 
ated by the oscillations will be the matter sector of the theory. If radiation 
domination is reached - and our model is predisposed to do so ^H] because 
radiation is the unique power law solution for which our simple source van- 
ishes (-R = 0) - the energy deposited into matter will sustain the radiation 
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domination. 

Each mode will contribute a very small amount of energy to the process 
but since there exist a huge amount of modes the whole process can be very 
efficient. Moreover, this reheating mechanism is very natural because the 
interaction of the modes with any matter field is of the universal gravita- 
tional strength. In contradistinction, scalar-driven theories involve different 
couplings of the infiaton to different matter fields. In one sentence, the scalar 
$ via its coherent oscillations can naturally reheat and lead to a hot thermal 
universe. 

The simple phenomenological model used in this paper has late time 
evolution problems which can be addressed by modifying the ansatz for the 
gravitationally induced source flH]. This change, however, only affects the 
late time evolution and does not disturb the results of this paper concerning 
primordial scalar perturbations. 
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